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Abstract 

We study elementary particle reactions that result from the interaction of an atomic system 
with a very intense laser wave of circular polarization. As a specific example, we calculate the 
rate for the laser-driven reaction e"'"e~ where the electron and positron originate from a 

positronium atom or, alternatively, from a nonrelativistic e'^e~ plasma. We distinguish accordingly 
between the coherent and incoherent channels of the process. Apart from numerical calculations, 
we derive by analytical means compact formulas for the corresponding reaction rates. The rate 
for the coherent channel in a laser field of circular polarization is shown to be damped because of 
the destructive interference of the partial waves that constitute the positronium ground-state wave 
packet. Conditions for the observation of the process via the dominant incoherent channel in a 
circularly polarized field are pointed out. 

PACS numbers: 13.66.De, 41.75.Jv, 36.10.Dr 
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I. INTRODUCTION 



The interaction of electrons and atoms with laser radiation is intensively and successfully 
being studied for many years now. However, due to a rapid technological progress, the high- 
power laser systems available today can generate peak intensities up to 10^^ W/cm^ in the 
range of near-optical infrared frequencies jH , and a further increase can be expected within 
the next few years Consequently, the ponderomotive energy of an e~ (or e"*") inside 
such a laser wave is of the order of 1 GeV, which is far beyond the typical energetic range 
of atomic physics but rather reaches the energy scale characteristic for elementary particle 
physics. If electrons or an e~ and collide at such high energies, then particle reactions 
like heavy lepton-pair creation or hadron production can occur. This indicates that there 
might be a way to merge laser physics with high-energy physics ISLyi. Similar efforts are 
being undertaken with respect to laser physics and nuclear physics j^, 0, 0, 13| • 

The high-energy process of photon-induced e^e~ pair creation by a projectile particle col- 
liding with an intense laser beam has already been investigated before, both experimentally 
and theoretically l^, Hi- An essential ingredient to these studies is the ultrarelativistic 



energy of the incoming particle. In its rest frame, the doppler-shifted laser frequency and 
field strength are considerably enhanced. As a consequence, the projectile actually faces an 
x-ray beam of near-critical intensity. 

Instead, in the present paper we study a situation where elementary particle reactions 
arise from the interaction of a strong laser field with a nonrelativistic atomic system. To this 
end, we suppose that a positronium (Ps) atom is brought into an intense laser wave. We 
note that the lifetimes of ortho-Ps (~ 10~^ sec) and para-Ps (~ 10~^° sec) are much longer 
than the typical duration of a strong laser pulse. Due to the equal masses of its constituents, 
the dynamical response of the positronium to the electromagnetic forces exerted by the laser 
field is rather unique 2|: The laser's linearly polarized electric field leads to an antiparallel 
oscillatory motion of the particles in the transverse direction, while the magnetic Lorentz 
force causes an identical ponderomotive drift motion along the laser pro pag ation direction. 
This leads to periodic e'^e~ (re) collisions (see, in particular. Fig. 1 in Ref. |l4|). If the energy 
of the relative e^e~ motion is large enough, then in these coherent collisions particle 
reactions can occur. Thus, we shall study high-energy processes induced by e~^e~ annihilation 
resulting from a laser-driven Ps atom. Considering the case of a circularly polarized laser 
field we will find as a main result, however, that the various partial waves that constitute 
the Ps ground-state interfere destructively, which causes a heavy suppression of the coherent 
reaction rate. This quantum effect can be related to the classical trajectories of the colliding 
particles in the laser field. Furthermore, when the characteristic size of these trajectories 
(or the size of the spreading particle wave packets) exceeds the interatomic distance, then 
collisions between particles originating from different Ps atoms will come into play, which 
opens the incoherent channel of the process. Surprisingly it turns out that, in a circularly 
polarized laser field, the incoherent channel is dominant as the interference in the coherent 
channel is destructive. In order to study the incoherent process we replace the Ps atom by a 
nonrelativistic e^e~ plasma. In this situation, exclusively incoherent e~^e~ collisions occur. 

It should be stressed that in the described setup the e~^e~ collision energy is basically 
determined by the kinetic energy ~ mc^^ contained in the transversal motion of the par- 
ticles, which is considerably smaller than the ponderomotive energy ~ mc^^^ mentioned 
above. Here, mc^ is the electron rest energy and ^ = ea/mc^ denotes the so-called laser 
intensity parameter with the electron charge — e and the laser's vector potential a. For the 
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highest intensities attainable at present ^ is of order 10^. In this respect, the underlying 
laser acceleration of the particles is considerably different from the usual laser acceleration 
techniques, since the latter try to extract the ponderomotive energy gain along the laser 
propagation direction. Nevertheless, the energetic thresholds for muon or pion production 
might be within reach. We further notice, that the principal difficulties of laser acceleration 
implied by the Lawson- Woodward theorem (see, e.g., Ref. 0|) are completely absent here 
since the e+ and e~ collide inside the laser wave. 

Against this background, we consider the specific process e^e^ A*^/^^ which is one 
of the most fundamental in high-energy physics. Its cross section sets the scale for all 
e'^e~ annihilation cross sections [l^. For example, at high energies one has CTg+e- ^hadrons ~ 
4cre+e-^/x+/x-5 where o"e+e-^hadrons denotes the total cross section for the production of any 
number of strongly interacting particles The threshold energy for the reaction e^e" 
fi~^fi~ amounts to 2Mc^ in the field-free case, where M denotes the muon mass. According 
to the above, a naive estimate thus suggests that a laser intensity corresponding to ,^ ^ 
M/m ~ 200 is required to produce a muon pair in a laser-driven e~^e~ collision. This value 
is reached, e.g., for a linearly polarized laser beam of 3.8 x 10^^ W/cm^ intensity and 1 eV 
photon energy. 

To the best of our knowledge, the process e~^e^ ^^fi~ in a laser field has not been 
considered before. The most closely related article treats the laser- assisted Bhabha scattering 



e ' e — >• e ' e 



17| . In Ref. the low-intensity case (i.e., ^ -C 1) is analysed in detail with 
the emphasis lying on the resonances that can occur in the scattering cross-section due to 
the interaction of the leptons with the background laser field. We will come back to this 
point later. Another similar process, that has found the interest of several authors, is the 
M0ller scattering e~e~ e~e~ in a laser field (see jisl . lol 2ol 21 1 and references therein). 



The paper is organized as follows. In Sec. II we develop a formalism that allows us to 
calculate the rate for the reaction Ps in a strong laser field. Our treatment will be 

based upon the Volkov solutions to the Dirac equation. Afterwards we analyse in detail the 
reaction kinematics. Here we show in particular that the minimal laser intensity parameter 
required is indeed given by ^min = M/m [cf. Eq. flHUjl ]. Further, the kinematical analysis 
will help us to derive a compact formula that gives an approximation to the total reaction 
rate and displays its main dependences [cf. Eq. (jHH)]- In Sec. Ill we present our (numerical) 
results on the total and differential production rates and compare them with the known 
cross section for the field-free process e~^e~ —>■ fi~^fJ'~ ■ Furthermore, we briefly consider the 
related process of muon pair production by a superstrong laser wave interacting with a 
nonrelativistic e~^e~ plasma. In this situation the interference effect does not play a role. 
We finish with a conclusion. 



II. THEORETICAL FRAMEWORK 

A. Transition amplitude and reaction rate 

We calculate the rate for positronium decay into muons in a strong laser field, i.e., the 
rate for the laser-driven process Ps fi^fi~ . We assume a photon energy of about 1 eV and 
a laser intensity parameter of order M/ m ~ 200 or larger ^2] • For mathematical simplicity, 
the laser field is taken to be a monochromatic, plane wave of circular polarization with the 
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classical four-potential [2 



A^{x) = tti cos{kx) + 02 sin(A;x) 



As usual, is assumed to be adiabatically switched on and off in the remote past and the 
distant future, respectively. In Eq. (jT)), = 0,0,1) is the wave four-vector and 2 
are constant four-vectors chosen as ai = (0, a, 0, 0) and 02 = (0, 0, a, 0) with a denoting the 
amplitude of the vector-potential. From now on we use relativistic units {h = c = 1), except 
where otherwise stated. We notice that in the circularly polarized laser field (HJ the e"*" and 
e~ are permanently colliding since, according to the classical equations of motion, they are 
co-rotating in the polarization plane. 

The Ps atom is assumed to be initially at rest and in its ground state. In a usual field 
theoretic formalism jl^, this bound initial state can be expressed as a superposition of 
products of free states ipp^ for the electron and positron with definite momenta = ±p. 
The superposition is weighted by the probability amplitude for finding a particular 

value of p. Note that this amplitude is just the Compton profile of the Ps ground state (i.e., 
the Fourier transform of its wave function) and p can be viewed as the relative momentum 
of the electron-positron two-body system (i.e., as the momentum of an effective particle of 
reduced mass m/2). When submitted to the strong laser field (^ ^ 200) the Ps atom will 
instantaneously be ionized, and the dynamics of the ionized e~ and e"*" will be governed by 
the laser field, which predominates over the influence of the Coulomb interaction between 
the particles. Therefore, in the spirit of the strong- field approximation theories jl^], we 
may replace the free leptonic states ipp^ by laser-dressed Volkov states 2^ 2^ . Within this 
framework, the amplitude for the laser-driven process Ps — > yU^/i~ can be written as 



VJ (27r)= 



$(p)5e4 



(2) 



with a normalization volume V and 



-iaf / d^x / d^y^p^,,^{x)Y^p_,s_{x)D^,{x-y)^p_,sAyh'''^P+,sAy) (3) 



being the amplitude for the process e+e — * fi^fi in a laser wave (cf. Fig. 1). In Eq. (jS)), af 
denotes the finestructure constant, 



D^uix -y) = j 



d^q ei-^''-y) 

(27r)4 g2 
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(4) 



is the free photon pr opagato r |27l l28l |. and the laser-dressed states for the electron and 
positron are given by |25l . |26[ 



^p±,s± (3^) 



^0 * 2(A;p±)j''"±'^±'' 




(5) 



with 



= ±{g^x) + sin(fcx) - cosikx) 



{kp±) 



{kp±) 
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In Eq. (0), p± are the initial free four- momenta of the electron and positron (out side the 
laser field), s± denote the particle spin states, the Mp±,s± are free Dirac spinors |29l|, and 



(6) 



are the effective four- momenta of the particles in the laser field ^I]. Note that Eq. © implies 
= and thus + = 0, where the label _L denotes the momentum component that 
is perpendicular to the laser propagation direction. The corresponding effective mass reads 
ml = = (1 + S,'^)m? with the dimensionless laser intensity parameter 



ea 
m 



(7) 



Like free states, the Volkov states in Eq. (0) are normalized to a (5-function in p± space 
26l 3^ . Analogous expressions hold for the Volkov states p^^s±i the free momenta P±, the 
spin states S±, the effective momenta the effective mass M, = M(l + S2)V2^ and the 
intensity parameter S = ea/M of the muons. Note that the amplitude (jHl) fully accounts 
for the interaction of the leptons with the laser field, while their interaction with the QED 
vacuum is taken into account to lowest order. Similar approaches have been used for the 



theoretical description of laser-assisted e+e" 113 and e-e" [Ulli,!!! I21l scattering. 

By the standard procedure of using the generating function of the Bessel functions (3ll |. 
one can perform the space-time integrations in Eq. (jSJ to get 



X ^{q+ + q 



m 



M 



/popo 

q - nk) 6{Q+ + -q- Nk) 
with the electronic spinor-matrix product 

2{kp^){kp.)) " 



(8) 



M''{p+,p-\n) 



r 



+ 



+ 



2{kp+) 2{kp^) J 
2{kp^) 2{kp_) ) 



(9) 



and a corresponding expression A4^{P^, P_\N) for the muons. The coefficients in Eq. Q 
are given by 



b: = 




a)e 


inipo 




bt = 


1 

2 




a)e~ 




b- = 


1 

2i 


Jn-1 


(a)e~ 





with a 



al + q;2, v^o = arccos(a;i/a) = arcsin(a2/a;), and 



a A 



{kp- 



{kp+) 



(10) 



(11) 
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for j = 1,2. As is expressed by the energy- momentum conserving 5- function at the first 
vertex, the integer number n in Eq. (jH} counts the laser photons that are emitted (if n > 0) 
or absorbed (if n < 0) by the electron and positron. Similarly, N is the number of laser 
photons emitted (if < 0) or absorbed (if > 0) by the muons. Denoting the total number 
of absorbed laser photons by r := A^ — n and integrating over the virtual photon momentum 
yields 



m M 

Se+e-^^.+f.- = -i(27r)V^==^==^A<^(p+,p„|n)A^^(P+,P_|n + r) 



^ 5iq+ + q--Q+-Q- + rk) 
(g+ + g_ — nk)'^ 

In general, the denominator (g+ + g_ — nk)'^ in Eq. ()12|) could become zero. By way of a 
renormalization procedure, such mathematical singularities can be transformed into physical 
resonances that appear in the production process 17, [H, [l^ 23|- easily see. 



however, that in the present situation, due to the large value of the laser intensity parameter 
and the nonrelativistic electron and positron momenta p±, one is always far off resonance 
js^l- Namely, on the one hand we have 

(g+ + Q- — nkY = 2ml + '^(q+q ) ~ 2n{kp^) — 2n{kpJ) ^ 4m^ — Antum 

which becomes zero for 

n.es ^ e- - 10^°. (13) 

On the other hand, the Bessel functions J„(a), that enter the production amplitude through 
the coefficients in Eq. (fTUj) . practically vanish unless a ^ n. Since g+ ~ q'^ and m^,q± ^ Qz, 
the argument approximately equals 

\q \\ . ^TTi\q \\qz 



uj go- Qz ' ml 

where we have dropped the particle labels ±. Now, \q^\ = \p±\ ~ wiaf and Qz ~ m^'^/2. 
Hence, 

a ~ Of^ 10^ (14) 

which, according to Eq. (|13p. is orders of magnitude smaller than would be required for a 
resonance to occur. 

The above argument can be further exploited. From Eq. p4j) we know that the main 
contribution to the production amplitude comes from photon numbers n with |?t,| ^ 10^. 
But for those numbers we have to a very good approximation 

(<?+ + Q- - ^ (g+ + qSf 
which, thus, can be pulled out of the sum in Eq. (fT^: 

m M 1 

5...-.,.,- - -i(2vr)V-^^^— ^EA^^XP.,P-l-)A^.(n^ 

X 5{q+ + q_-Q+-Q_+rk). (15) 
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The summation over n can now be performed analytically by virtue of Graf's addition 
theorem l3lll with the result 



J2M^{p+,p.\n)M^{P+,P.\n + r) = J^u^U^ + K+iu^'V^ + v'^U^) + K; {u'^W^ + w'^U^) 

n 

+ Ltv'^V^ + Mriv'^W^ + w^V^) + L-w'^W^. (16) 

Here we have used the abbreviations 



u 



u 



u 



r- 




2{kp+){kp-] 






2{kp. 


f) 2(A;p_) 








f) 2(A;p^) 



V— T^- 



u 



P— JS- 



u 



V— I*- 



(17) 



and similarly U^,V^, and for the muons. The coefficients in Eq. (fTHIl read 









i 

~ 2 






k; 


i 

~ 2 








1 

~ 4 


2Jr{Sy ± Jr^2i5y~^ ± Jr+2i5) 


Mr 


i 

~ 4 







,r+2 



(18) 



with 



7 = /? — ae 



i(¥'o-'?o) 



and e = ^e*^° 





where [3 and ?7o are the muonic quantities that correspond to a and ipo. We will see later 
that 7 ~ /3 since a ^ j3 for the typical parameters. An insignificant overall phase factor of 
e'^''° can be dropped in Eq. (fTHj) . 

Now we come back to the reaction Ps In order to obtain the corresponding 

amplitude we have, according to Eq. (0), to multiply Eq. (fT3j) by the Compton profile $(p) 
of the positronium ground state and integrate over the relative momentum p. It turns out 
that this integration is a very difficult task that can only be done in an approximate way: 

First, within the momentum range given by $(p) the electronic spinor- matrix products 
in Eq. ()17|) are practically constant (on the 1% level since \p\/m ~ af) and can therefore be 
pulled out of the integration. The same holds for the kinematic factors p^. ~ m, (g+ + g_)^ ~ 
4m^, and the energy- momentum conserving (5-function js^]- Hence, we are left with integrals 
of the form 



1 



(19) 



where x = arctan[asin(93o — ''?o)/(q^cos(9?o — "^^o) ~ P)] such that exp(ix) = 7/5. The highly 
oscillating factor exp(irx) leads to a very small value of J^. The oscillatory damping of the 
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amplitude is due to a destructive interference of the various partial waves within the Ps wave 
packet. Classically, this interference effect can be related to the extended motion of the e"*" 
and e~ in the polarization plane of the laser. Therefore, the mean impact parameter of the 
e'^e~ collisions is much larger than the initial Ps size and the resulting /i"*"//" production 
amplitude is suppressed. The laborious evaluation of the integral (fT^ is performed in the 
appendix. The result is 



Jr ^ ( ^Yp'^'MP). (20) 



with the Ps radius = 2/ aim. The damping factor can also be written as (tu/maf^)^ = 
(vrao/A^)^ ~ 10~^^. Note that 2A^ gives the average impact parameter of the e"'"e~ collisions 
since, according to their classical trajectories, the particles co-rotate in the polarization 
plane on opposite sides of a circle of radius A,^. However, the classical picture suggests 
that the process probability is proportional to (ao/A^)^. Instead, this factor is contained in 
the process amplitude such that the probability scales as (ao/A^)^. This indicates that the 
damping factor is truely of quantum mechanical origin. 
The square of the amplitude reads 



J2M''{p+,P-\n)Mf,{P+, P-\n + r) 

n 

X 6{q+ + q--Q+-Q-+rk)VT, (21) 



where J2n indicates the sum over n in Eq. (fTI)|) averaged over the Ps ground state, as described 
above, p± {q±) are to be understood as some typical values of the electron and positron 
(effective) momenta, and the factors of volume V and time T come, as usual, from the square 
of the 5-function. Note that the energy-momentum conserving 5-function, in particular, 
implies + = 0. From Eq. (PT|) we get the total reaction rate by averaging over 
the initial spin states, summing over the final spin states, and integrating over the final 
momenta: 

_ 1 rd^P^ rd'P^ 1 

Rps^,-^,- - ^ y y 4 ^2^^ \sps^,^,- 1 . (22) 

In the next but one subsection we derive a compact analytical formula that gives an estimate 
for the muon production rate (j22p . But before, we analyse in some detail the reaction 
kinematics. 



B. Kinematical considerations 

In the following we provide estimates for the minimal (rmin) and the typical (f) photon 
numbers that are net-absorbed during the production process. From the latter we also find 
the typical momenta of the created muons. 

A lower bound on r can be derived from the equation 

(g+ + g_ + rkf = 2(g+ + g_ + rk) ■ Q± (23) 
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that follows from the energy-momentum conservation condition expressed by the 5 function 
in Eq. (j^H). Setting Qr = + + ru and qr = it ~^ 1z ~^ '^'^) this can be rewritten as 



cos^^r 



2QrQl - {Q'r - ql 



2gr|Q±| 

with the polar angle 6q^ = Z(fc, Q±). Demanding cos^ 9q^ < 1, we get 



(24) 



Or 

2 




AM: 



2 \ 1/2 



(25) 



Hence, it is required that 4M^ < Ql — q^ = {q+ + q- + rk^, i.e., the laser-dressed collision 
energy has to exceed twice the laser-dressed muon mass. Using (g++g_+rA;)^ ^ Aml+ircum, 
we find 



^ ~ ^min 



M2 - M2 



com 



(26) 



This means that, e.g., for uj = 1 eV at least 2 x 10^'^ photons have to be absorbed from the 
laser wave for muon production to take place from the initially low-energy e^e~ pair. This 
number is independent of the laser intensity. 

Assuming a symmetric situation with ^ Q*^ , Eq. ()26|) implies that the minimal muon 
energy is approximately given by 



(27) 



Hence, the muons are typically produced with highly relativistic momenta such that their 
dispersion relation approximately reads Q\ ~ |Q±|- Furthermore, they are emitted roughly 
along the laser propagation direction (note that Qz ^ Q± since q^ ~ m^^ -|- ra;, Qr ~ 
2m + m^^ -|- ruo so that q^ ~ Qr)- More precisely, by solving Eq. (j24j) for Q\ we find that 
the polar emission angle satisfies the relation 



cos6'q^ > — 



Qr 

qr 



1 - 



{Ql 



qlf 



AQIMI 



1/2 



AvfT 



1/2 



(28) 



For these reasons one can say that the muon kinematics is similar to that of the laser pho- 
tons. This "photon-like" nature of the muons results from the fact that they are essentially 
produced by a huge number of laser photons whose total energy, according to Eq. (jTTjl . 
exceeds the initial nonrelativistic energy of the e~ and e"*" by orders of magnitude. 

The typical number of absorbed laser photons can be estimated by exploiting the prop- 
erties of the Bessel function Jr(/5) in Eq. (j2(J|) . To this end, let us again assume a symmetric 
situation, which allows us to drop the particle indices ± in what follows. The energy- 
conservation condition then can approximately be written as IQq ~ 2go + tuj. Because of 
the photon-like muon momenta this can be expressed as 



2Q, 



'2QI 



2qJl + 



2ql , 



+ ruj, 



(29) 
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where ^ qo ^ Applying the momentum conservation condition 2Qz ~ 2qz + ru, 

we thus get 

Ql + ^ ql + ml 

Now, let ro = q^/uj and ^ = t/tq. Then, again by the momentum-conservation condition, 
^ (1 + i/2)q^. Hence, Eq. ^ implies 



1/2 



(31) 



where m* gj, was used. Thus, the argument of the Bessel functions in Eq. (j20| approxi- 
mately equals 

1 /2 

ea Q\ m Q\Qz Q\ m ( ru , . 

^ ^ 2- n n ^ ^^- n^i A!f2 ^ 2^— ^ 2^- ^ . 32 

^ Qo-Qz oj + UJ oj \m ) 

We note that, according to Eq. ()2fi|l . the expression under the square root on the right-hand 
side of Eq. flH^ is positive. By the properties of the Bessel functions [Slj], the typical number 
of absorbed laser photons is expected to be determined by the condition (3 k, r. This yields 



r ^f = 2i''-(l + Vl^) (33) 
with K = M/m^ ji^l- The corresponding typical muon momenta read 

a « |^2[i^2(l + v^r^)]. (34) 



Using the Eq. ^ between the effective and the free four-momenta and the relations ()45|1 
below, we find for the typical values of the muon momenta after the interaction with the 
laser field 

= Ox ~ K (i + v^r^) , 

Pz ^ Qz- ~ me (i + VT^) . (35) 

For example, for ^ = 250 and = 1 eV we have f ^ 10^^, Q± = P± ^ 2M, Qz ~ 620M, 
and Pz ~ 470M. From Eqs. ()33|1 and ()35|1 we see that the typical final energy of the muon 
pair satisfies the relation 2Pq ^ fu, which reflects the law of energy conservation after the 
laser has been switched off. Equation (|33p also implies that the minimal intensity parameter 
required for the process to have a significant probability (i.e., to be able to fulfill r ^ P{r)) 
amounts to 

^min = — (36) 
m 

which agrees with our earlier naive estimate. 



10 



We notice that the partial rate for muon production by the absorption of r = rmin 
photons is zero. Namely, according to Eqs. (jHT|) and for £ = rmin/To ~ 2M^/m^^^ the 
transverse muon momentum and with it the argument of the Bessel functions practically 
vanish. According to the above, the partial rate reaches a maximum at r ~ f, which at 
^ = ^min is twice as large as the minimal number of photons: f = 2r^i^. 

It is interesting to observe that the typical muon momenta in Eq. (j35|) can be interpreted 
by employing a classical simple man's model of the creation process. In the classical picture, 
the threshold value of the laser intensity (jHUj) corresponds to the situation when, in the 
center-of-mass frame of the e~^e~ system, the kinetic energy is large enough to create muons 
at rest. I.e., denoting the laser phase by r = Lj(t — z), we have 

P'±iro) = P'Aro) = (37) 

at the creation phase tq, where the prime indicates the center-of-mass frame. The classical 
equation of motion for a muon in a laser field with the initial condition (j37p has the solution 

Pi(r) = eA(r) - eA(ro) , 

P'Ar) = ^Kr)-A(ro)]^ (38) 

Consequently, after the interaction with the laser field the muon momenta equal 

P[ = m^, 

P' = (39) 

Due to the e+e~ longitudinal drift motion in the laser field, the relative velocity between the 
center-of-mass frame and the lab frame amounts to frei = <iz/<lo = ^^/(2 + ■C^) [see Eq. (jH))]. 
The Lorentz transformation to the lab frame thus yields 



Pj_ = mi, 



2 



^'-^TiTF + ij^Vi+e (40) 

The latter coincides with the typical muon momenta at the threshold ^ = ^min given by the 
quantum theory: k, M"^ /m and P^ k, M [see Eq. ()35p ]. The typical number of absorbed 
photons is determined by the muon final energy: fuo = 2{Pq — m) ^ 2M^/m, which is in 
agreement with Eq. (jHEI)- 

Our simple man's model can also explain a peculiarity in the angular distribution of the 
muons (see Fig. 4 in Sec. III. A). Since Pz ^ P±, the muons move in a narrow cone with 
the axis parallel to the laser propagation direction, but at very small angles the angular 
spectrum has a dip. The dark region in the angular distribution occurs because the muons, 
although having been created with zero transverse momentum in the laser field, acquire a 
nonvanishing transverse momentum after switching off the laser field [see Eqs. ()37|) and 



C. An approximative formula for the total rate 

Equation (j^^ for the total rate of the reaction Ps f^ f^^, although looking rather 
innocent, is actually quite involved and can be evaluated only numerically. Therefore it is 
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desirable to find, by analytical means, an approximation to Eq. ()22p that displays its main 
physical content. To this end, we consider the contribution to the total rate stemming from 
the first term on the right-hand side of Eq. ()16|) . From our numerical calculations we learn 
that this term gives by far the main contribution (~ 90 %). Thus, we need to calculate 

X J2W)?(^'^' E \u^U,\'6{q^ + q^-Q+-Q^+rk). (41) 

r s±,S± 



With the help of the (5-function and the relations |2 

^ = ^, d'Q^ = \Q_\QUQUcoseQ_ 
we can integrate over d^Q^ and dQ^ to find 



(42) 



where the relations ~ m, (g+ + g_)^ ^ 4m^, and ~ have been used. The spin 
sum in Eq. fjl^ can be converted in the usual way into a product of two traces: 

%iU,uU ■= E I'^^^t^f 

s±,S± 

^ j ~ 2{kp+){kp_) ) 2m y ~ 2{kp+){kp^) J 2m J 

2{kp+){kp_)) 2M 2{kp+){kp_)) 2M ]■ ^ ^ 

We want to find some typical value of Tuu,uU- The standard trace technology yields 

[{P-P-) + iP-P+) + {P+P+) + ip+P-) - 2{P+P-) - 2(p+p-)] 

Here we have used the (remarkable) relations 

um ^ {kp+) ^ (kp^) ^ {kP+) ^ (kP^) (45) 
that hold to a good approximation since 

(kP) = (kQ) ^ ^ ^'^i^ ^ = 

by Eq. (jHUjl . With (p+P-) ~ m^, (p±P±) ~ rnP^, and (P+P^) ^ + 2P| the expression in 
Eq. becomes 

W-^(Po-meY + 8e^§ (47) 
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with Pq and P± denoting some characteristic values of the muonic energies and transversal 
momenta that, by Eq. amount to Pq ~ 2m^^ and P± ~ 2m^ at ^ ^ ^min. This leads 
to the desired typical value of 



.4 m'^ 



Tuu,uu ~ see"—, (48) 



which can be pulled out of the integration in Eq. ()42|) . We proceed by performing the further 
approximations 



Jd<PQ_^2n, |cicosV^%, (3'/'^f'/\ J2iMm)?^^, (49) 
where, according to Eq. ()28j) . the maximum polar emission angle is given by 

« (50) 

Putting all pieces together, we arrive at the handy formula 

~ 32 ( u? \^ (^me\^"' 1 ,,,, 



where r^, denotes the classical electron radius [35[- Equation (fFT|) is the desired analytical 
estimate for the rate of laser-driven Ps decay into muons. 
We notice that Eq. (jSip can also be represented in the form 

Here a = {9/8){a^/M^) stands for the process cross section, which for ^ ^ ^^i^ becomes 



(53) 



Recalling the simple man's model and using the electron energy in the center-of-mass frame 
p'o ~ ^6 can infer that Eq. ()52|) is based on the process cross section 



CT ~ ^, (54) 

T 

with 7 = Pq/tti ^ e being the electron gamma-factor in the center-of-mass frame. Equation 
(j5l|l is in accordance with the known field-free cross section for muon production in e^e~ 
collisions [see Eq. 

If the interaction volume contains positronium atoms, then the rate will increase 
correspondingly: 

<W=^^Ps^M+.- (55) 

Here we have assumed that each Ps atom independently creates a muon pair, i.e. there is 
no interference between electrons (positrons) stemming from different Ps atoms. The latter 
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is the case when the electron wave-packets from different atoms do not overlap, i.e. when 
-C where n is the Ps density. Otherwise, when the spatial extension of the electron 

wave-packet is large, then the gas of Ps atoms transforms into an e~^e~ plasma. Against this 
background, in the following we will denote the rate in Eq. (j55j) resp. (|22p as the coherent 
rate for muon production since the colliding e"*" and e~ originate from one and the same Ps 
atom. In contrast to that, the rate for muon creation from an e~^e~ plasma, that might have 
been formed from an initial Ps gas, will be refered to as incoherent rate. In this situation, 
electrons and positrons from different Ps atoms can collide which gives a total number of iV^ 
incoherent collisions. In the next section we will present our results both on the coherent 
and the incoherent channel of muon production. 

III. RESULTS AND DISCUSSION 

A. Muon pair creation by a laser-driven Ps atom: The coherent process 

Based on Eq. ()22|) . we have numerically calculated the coherent rate for /i+^u^ pair creation 
from a single Ps atom submitted to a strong laser field of circular polarization. The laser 
frequency has been taken to be = 1 eV, throughout. For the laser intensity parameter 
we have chosen the three different values ^ = 250, 500, and 1000. The corresponding 
laser intensities amount to 1.1 x 10^^ W/cm^, 4.5 x 10^^ W/cm^, and 1.8 x 10^^ W/cm^, 
respectively. 

In Fig. 2, the dependence of the total muon creation rate on the laser intensity parameter 
is shown. For the intensity parameters under consideration we find production rates of 
1.0 X 10-15 sec-i = 250), 1.6 x lO'^^ sec'^ = 500), and 1.1 x 10"^^ sec-^ = 1000). 
The analytical approximation (j3T|) overestimates these numbers, but is still in rather good 
agreement with them (cf. Fig. 2). The reason for the overestimation, in particular for 
^ ~ iCmin, are the rather large values of Pq and P± used in Eq. (jTTj) . A total production rate 
of 10~i5 sec~i means that in a finite laser pulse of femtosecond duration the probability 
to create a muon pair from a single Ps atom is of order 10"'^'^. We notice that a typical 
laser focal volume will contain only one Ps atom on average since the highest positronium 
densities achievable at present are on the order of 10*^ cm~^ j36|. However, proposals to reach 
a Ps density of 10^^ cm~^ js^l or even a Ps Bose-Einstein condensate of 10^^ cm~'^ |38| are 
being considered. The above numbers of created muons seem too small to be experimentally 
accessible. Clearly, the main reason for the smallness of the coherent reaction rate lies in 
the damping factor (ao/^O^ ~ 10~^^ [cf. Eq. (PU|) ]. The latter results from the destructive 
interference of the partial waves constituting the Ps wave packet in the laser field, which 
makes the recollision of an e"*" and e~ from the same Ps atom highly unlikely. From a 
simplified, classical point of view, the reason for the damping lies in the large collisional 
impact parameter of order A£ that is due to the equal handed rotations of the particles in 
the laser polarization plane [39|. 

In Fig. 3, the partial production rates with respect to the number r of absorbed laser 
photons are shown (i.e., the contributions to the total rate stemming from a net absorption 
of r laser photons in the production process). The photon number is given in units of 
ro ^ m^V2cu, which amounts to 1.6x10^° = 250), 6.4x10^° = 500), and2.6xl0ii = 
1000), respectively. The shape of the curves in Fig. 3 can be understood with the help of the 
kinematical analysis in Sec. II. B. First, according to Eq. (j26|) . the minimal number of laser 
photons required from kinematical constraints amounts to = 2.2 x 10^°, independent 
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of the laser intensity. If we express this number with respect to the respective values of ro, 
then we get r^in/ro = 1.4 = 250), 0.3 = 500), and 0.1 = 1000). The partial reaction 
rate for r = rmin is always zero, as can be seen in Fig. 3. Further, in agreement with Eq. ()33p. 
the curves exhibit maxima at r f; for ^ = 250 the maximum is located at t/tq ^ 6, while 
for ^ = 500 and 1000 we have f/ro ~ 8. This feature reflects the mathematical properties 
of the Bessel functions and is correctly predicted by our simple man's model. 

Figure 4 shows the angular distribution dR/dOq for one of the produced muons. For 
symmetry reasons, the spectra for the muon /i~ and the antimuon /i^ are identical. The 
differential rate is expressed with respect to the polar angle 6q = l{k,Q) of the effective 
muon momentum. The value of the laser intensity parameter is chosen to be ^ = 250. 
One can see that the muon is emitted into a very narrow angular range starting from 
1.3 X 10~^rad ~ 0.07° and extending to 4.7 x 10~'^rad ~ 0.27°, which is in agreement 
with Eq. (j50p . In other words, as has already been mentioned before, the muon moves 
practically parallel to the propagation direction of the laser beam. The occurence of a 
minimal emission angle arises from the fact that, according to Eq. (j32j) . the argument of the 
Bessel function is proportional to the transverse momentum component Q±, which itself is 
proportional to sin^Q ^ 6q. Thus, the emission angle cannot be too small because otherwise 
the Bessel function will vanish. Since P± = Q±, a dark angular region also exists in the 
angular spectrum dR/dOp with respect to the muon momentum outside the laser beam. An 
alternative explanation of this phenomenon in terms of a simple man's model has been given 
at the end of Sec. II. B. 



B. Comparison with the field-free process e^e ^ 

In this subsection we want to draw a comparison between the coherent muon creation 
from a laser-driven Ps atom and the corresponding field-free process e^e~ ^^fi~. In the 
high-energy limit, the cross section for this reaction reads [1| 



Air m? 2 

afrec = (56) 

o S 

where ^/s ^ 2M denotes the collision energy. In the case with laser field, the square root 
of the quantity 

(g+ + g_ + fkf ^ Ami + ^ruom ^ 20m^^^ (57) 

can be regarded as some average collision energy. For the ^-values considered in this paper, 
the laser-dressed collision energy is thus about 1 GeV. Hence, the reference cross section 
in Eq. ()56|) to compare with should be taken at ~ 1 GeV, where its value is about 
100 nbarn. To transform this cross section into a reaction rate, we have to multipy by the 
incident particle flux. In collider experiments, instead of the incoming flux, the luminosity 
is more commonly used. When a beam of iV+ positrons collides at high energy with a beam 
of electrons, then the luminosity is given by 

where U is the circumference of the collider ring and A is the beam cross sectional area 
at the collision point. To make the comparison with a single laser-driven Ps atom, we 
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use N± = 1 along with the typical values U ~ 10^ m and A ^ 10~^ cm^. Note that 
the corresponding mean impact parameter ~ ^/A of the field-free collision is of the same 
order of magnitude as the electron-positron spatial separation ~ in the laser wave. The 
resulting luminosity L ~ 10^^ cm~^sec~^ leads to a muon creation rate of 10"^'' sec~^. This 
number is considerably smaller than the production rates we found in Sec. III. A. However, 
in a real collider experiment one has bunches of N± ~ 10^° particles leading to much higher 
luminosities and reaction rates, of course. 



C. Muon production from laser-plasma interaction: The incoherent process 

We have seen in Sec. III. A that the large electron-positron wave-packet size during the 
motion in the laser field suppresses the coherent reaction rate dramatically. To reduce this 
size and achieve e~^e~ collisions at microscopic impact parameters, one can think of em- 
ploying different, more complicated field configurations 'IS^. Otherwise, the muon creation 
from Ps atoms will be dominated by the incoherent production channel introduced at the 
end of Sec. II. C. The latter coincides with the process of muon creation, when a low-energy 
e~^e~ plasma interacts with a strong laser beam. We have redone our calculation for this 
situation, i.e., for a free, initially nonrelativistic and e~ that collide in a strong laser 
field and create a pair by annihilation. Our results on the partial production rates 

for reasonable experimental parameters (see below) are shown in Fig. 5. The shape of the 
curves is similar to those in Fig. 3. 

From Eqs. (jSip and ()20|) one can infer that the total rate for muon creation from laser- 
plasma interaction approximately reads 

9 m^lN^N. 
16vr \4nt ^ ^ 

where N± denotes the number of electrons and positrons in the interaction volume V^nt, 
which is given by the laser focal volume. From Eq. ()59|) we can estimate the total number 
of produced muons A^^ = RrNg during the interaction with A^,, laser shots, each single shot 
having a pulse duration of r. When plugging in some reasonable numbers: r = 100 fs, 
V^nt ~ (lOA)^ ~ 10~^ cm'^, and assuming that the presently achievable number of positrons 
A^+ ^ 10^ ^ A^_ 40] can be compressed into the interaction volume or, alternatively. 



is created via a newly emerging laser-based technique |41[, then we get A^ ~ 1 muon 
production event during A^^ = 10^° shots. This number, being based on rather optimistic 
experimental parameters (especially concerning the positron compression), indicates that 
the realisation of the incoherent muon creation process might be not inhibitory difficult, but 
still it will be very hard with modern experimental techniques. 



IV. CONCLUSION 

In this paper we have studied ^'^ production by e^e~ annihilation from a laser-driven 
Ps atom. To this end, a calculational framework has been developed where the initial bound 
state is described as a superposition of Volkov states weighted by the Compton profile of 
the Ps ground state. By virtue of the interaction with the QED vacuum, which is treated in 
the first order of perturbation theory, this initial state can decay into a laser-dressed fi~^fi~ 
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pair. Also, the related process of muon creation by the interaction of a strong laser field 
with a low-energy e^e~ plasma has been examined. 

We have found that the minimal laser intensity required for the process to occur corre- 
sponds to an intensity parameter of ^ ~ M/m ^ 200. In the case of a near-optical laser wave 
of circular polarization, e.g., this value is reached for an intensity of 7 x 10^^ W/cm^. This 
means that, starting from a nonrelativistic Ps atom or e^e~ plasma, fundamental particle 
reactions can be ignited by a superstrong laser field of an intensity that is just one order of 
magnitude larger than the highest values available today. 

However, in the Ps case, the total production rate resulting from the coherent recoUisions 
is extremely small and amounts to about 10^^^ per second only. The strong suppression is 
caused by a destructive interference of the different partial waves constituting the bound 
initial state in the superintense laser field. This phenomenon is also expressed by a compact 
formula for the total rate that we derived by analytical means. As a consequence, in the 
considered setup the production rate will be dominated by the muon creation via incoherent 
e~^e~ scattering, for which the system of Ps atoms has no advantage compared to an e^e~ 
plasma. For the incoherent collisions, the interference plays no role and the resulting muon 
creation rate is significantly larger than the corresponding rate from the coherent production 
channel. Nevertheless, very demanding experimental conditions are required in order to 
achieve observable muon yields. 

Finally, we note that similar reaction rates can be expected for laser-driven tt+tt^ pro- 
duction from e^e~ or A more promising alternative within circularly polarized field 
configurations might be the process e~^e~ — e^e~ + e^e~, the field-free cross section of 
which being several orders of magnitude larger than the one for creation in Eq. (j56p 

i3- 
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APPENDIX 

In this appendix we calculate the integral in Eq. (fT^ . Using cylindrical coordinates it reads 

1 /•+00 



(27r)W 



dpz / p±dpj 
Jo 



d(p^{p)Jr{S) e 



(Al) 



with 



X 



+ 13"^ — 2aj3 cos(v9o ~ Vo) 

uj p]_ + rn? ' 

asin((/)o - ??o) 1 



1/2 



arctan 



[acos(v9o -rjo) - 13]' 



[1 + {a,p^? + {a,pz?f 



and (3 ^ r being of order 10^^ for the laser parameters under consideration [see Eq. (jHS))]- 
Note that ip^ coincides with the azimuthal angle if. Going over to the variables cioP± ~^ 



aoPz 



and {(fo — rjo) Eq. ()A1|) reads 



o/o ^ ; dz I xdx 



dip ■ 



Jr{6)e 



irx 



1 + + Z 



2^2 ' 



(A2) 



For brevity, the overall factor 2/ {n^aQVy^^ will be dropped in what follows and only restored 
in the final result [see Eq. ()A14|) ]. Since a <^ (3 we can expand the functions 6 and x'- 



a 



5 = (3 — a cos (p + O \ — 



a 



a 



X 



-sinc^ + Ol^cosV). 



Note here that, since 0^/(3 ^ 1 [see Eq. ()14|l ]. we may drop terms of this order in the 
expansion of 5. Further, since the Bessel function Jr{5) is exponentially or oscillatorily 
damped unless |r — (5| ^ ^-i/s^ ^^j^g main contribution to the integral comes from the region 
with |acosv5| ^ r^^^ ~ (3^^'^. Accordingly, terms of order r(a cos y?//?)^ ^ may be 

neglected in the phase rx- Hence, we obtain 



dz 



10 Jo (1 + + 
Next we observe that Jr = 2 Re 1^ with 

3y 



dip Jr{(3 — a cos ip) e 



la sin ip 



dip Jr{(3 — a COS ip) e~ 



■lasin ip 



OO /-OO 

Jo (1 + a;^ + z'^Y Jo 
Jo [\ + x^ + z'^y Jo 



^ Jo ^ Jo (1 + a;2 + ^2)2 



dip Jr{P — a COS ip) e" 



lasin ip 



(A3) 
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In the last step we exploited the fact that the contributing ranges within the exponential 
and oscillatory regions of the Bessel function have a similar size. Let now 



x^x^ + z'^ + 
C = 2 I -2 ' P = C cos 



with e = af/2, such that a = ao( with ao = 2m^/u}. Performing the substitution of variables 

(x, 2;, if) (x, C, p) we obtain 



Ir = 2 / 

Jo 



°° x^dx 



X 



+ €-2)3 J^^.^ 1.2 x2 r>2 _ x^e-'^^l) 



vc^^ 



J,(/?-aop)e-'°»V^^ 



with Cmin = x/ + e 2. As mentioned above, the value of — aop) will be exponentially 
small unless p ^ po = /^^^'^/cto- Therefore, we can approximately write 



2J.(/3) 



°° x^ dx 



/■2 _ £ 



2(£-2-l) 
(x2+e-2)2 



(ip 







(A4) 



with pmax = min{C,Po}. 

In the following, the integral in Eq. ()A4|) will be evaluated in several steps. Let us first 
consider the integral over p. We show: 



/■Pmax p-iaoVC^-V 



TT 

2 



TT e 4 



-i«oC 



2 VC^ 



I C 



Note that, for the parameters of interest, we have 
aoPo ~ 10~2. In the first range ^ c^q "'^) we get 



-1 



(C«"o'). 

(ao ' < C < aoPo), (A5) 
(C > aoPo)- 

10"® and po ~ 10"^ such that 



C dp 



TT 







The second range we split into two cases. For q;q ^ <^ ^ < po the integrand is highly 
oscillating and the main contribution comes from the boundary terms: 



/(C) 



- i«oC 



+ / dy 



-iooC 



dt 



V2C Jo V 2 ./C^ 



-i«oC 
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where in the first and second steps the substitutions y = a/ ('^ — p'^ and t = ( — y were made. 
Note here that the contribution from the lower bound can be neglected since OqC ^ 1- 
Similarly, for po < C c^oPo we can write 



HO 



■ dy- 



-laoy 



+ 



we-Pi VC - iaopo 

Finally, in the third range (C ^ ctoPo) obtain 



dy 



-laoy 



/(C) ^/ fe 
Jo C 



C 



This shows Eq. ()A5|) . 

Now we continue with the integrations over dx d(. To this end, let us divide the integra- 
tion range into five regions: 



with 



I 


< X < Xi , 


Cmin < C ^ Cl 


II 


< X < Xi , 


Cl < C < C2 


III 


Xl < X < X2 , 


Cmin < C ^ C2 


IV 


< X < X2 , 


C>C2 


V 


X > X2 , 


C — Cmin 


X2 = 


«OPo T ^ 

~ -L , Cl 

e 


= -~io-«, 

ao 



C2 = "oPo ~ 10" 



Note that Ci and (2 coincide with the values of Cmin taken at positions xi and X2, respectively. 
Hence, it is easily seen, that this division covers the whole range of integration, i.e.: 

Taking into account Eq. ()A5|) . we find in the first range 

n(aoe)-i X^dx /-"o^ (dC 



x^eHt 



Jr{l3)TT^ I X* dx 

"'0 Jx 



, vc^ - x^e' le 

y/C^ — x'^ (C^ — + x'^e'^y 



-i)r 



where in the second step the transformations a^ex —>■ x, aoC ~^ C have been made. Intro- 
ducing the variable t = ("^ — x'^, this becomes 



TT 




2 


^/ 


TT 


1 


2 


age2 



dt 



x/t(t + x2e2)2 

(r^^^T^ + i) 



xe 



(A6) 



where in the second step formula 1.2.15.13 from Ref. |43|] was used. Since both integrals in 
Eq. ()A6|) are of order unity, we conclude that 



2 a^e2 ■ 



(A7) 
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Turning to the second range, we get similar as before 



The main contribution to the highly oscillating integral over comes from the lower bound- 
ary. In this way we find 

Jo i./T^{l - + x^e^y ^ ^ 

The remaining integral can be done analytically. Its value is ^e~^ such that we get 

^\ 3/2 g-i^ g-i 



ir'^MP){^) (A9) 



In the third region the integral reads 

= MP) 



e /""oPo 4 _.. , /-"o^o-^ y/t + x e"'* dt 



Jr{P)V2^-%e'^ r'^x^e-'^dx T 

"'1 "'0 



t(t + 2x) [t{t + 2x) + x2e2 



where first the same transformations as before were made {a^ex — x, aoC ~^ C) cind then 
the variable t = ( — x was introduced. Because of the highly oscillating integrand and the 
singularity, the main contribution to the t integration comes from the region around t = 0. 



Thus, we arrive at 



tto ^1 V2^{xey Jo Vi 



MP)^ (e--oPg_e-). (AlO) 



According to Eq. ()A5j) . in the forth range we have 



jy^) = 2J,(/3)po / 

JO 



2Jr(/3)eQ!oPo y x'^ dx J 



where we have substituted {aop^e ^) ^x ^ x and (aoPo) ~^ C- The highly oscillating 
integral over ( can be evaluated with the same methods as before to give 

JUV> « 2M,3) f (All) 

«o "'0 Wl - x^ [1 - x^ + x^e^)^ 
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The integral over x we have aheady encountered in Eq. ()A8|1 . Hence: 

^ ^r(/5)^^e-'"oP§. (A12) 

In the fifth range, after the substitution ex — > x, the integral can be written as 

dx /■°° e-^"o'^ rfC 



2 • 



With the definitions a = x/ a/1 + and 6 = x^(e^ + x^) /(I + x^)^, the ( integral can be cast 
into the form 



a VC' - a' -0? + by V2^b^Jo Vi V 2^0 v^fe' 



g-i«oa 



where the substitution t = ( — a was made and the oscillatory nature of the integrand 
exploited. This yields 

The main contribution to this integral stems from the lower boundary since for x ^ 1, where 
the oscillations fade away, the integrand is damped by the power x^^^/^. For this reason we 
obtain 

/ — e 4 e oPo 

/f) ^ MP)V2^ • (A13) 

In summary we have shown that 

i.e., the main contribution to the integral ()A4|) comes from region IV. Consequently, 

I ^2Reir = - 3,, ?/2 ^ ^ sin (A14) 

where we have restored the factor 2/ (rc^aQVY^'^ that was dropped from Eq. ()A2|) . 

We note that Eq. ()A14|) introduces a fast oscillating factor sin^(/3^/'^) to the differential re- 
action rate. Since (3 ~ 2m^'^ /uo these oscillations depend on the laser intensity and frequency 
[see also Eq. ()32|)]. In reality, however, neither of these parameters has a definite value in 
a short laser pulse but spreads over a certain range. In an experiment the oscillations are 
therefore averaged out, and we shall do the same: 
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FIGURE CAPTIONS 



Fig. 1. Feynman diagram for the process e+e in a laser field. The arrows are 

labelled by the free and effective momenta of the corresponding particle [cf. Eq. ©]. 

Fig. 2. Total rates for the laser-driven reaction Ps — >■ as a function of the intensity 

parameter ^ of the applied laser field. The black squares show the results of our numerical 
calculations based on Eq. (j221); the solid line shows the analytical estimate via Eq. (jSH). 

Fig. 3. Partial rates for the laser-driven reaction Ps — >■ /i"*"/^" as a function of the number 
of absorbed laser photons r for an intensity parameter of ^ = 250 (solid line), 500 (dashed 
line), and 1000 (dotted line). The latter two curves are enhanced by factors of 10^ and 
5 X 10^, respectively. 

Fig. 4. Angular spectrum with respect to the polar angle of the effective momentum for 
one of the produced muons at ^ = 250. 

Fig. 5. Partial rates for the reaction e"'"e~ fi'^fi~ from a laser-driven nonrelativistic 
plasma. It is assumed that N± = 10'' particles are contained in the interaction volume 
Vi^t = 10~^ cm^. The laser intensity parameter is ^ = 250 (solid line), 500 (dashed line), 
and 1000 (dotted line). The latter two curves are enhanced by factors of 2 x 10^ and 2 x 10^, 
respectively. 
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